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Abstract 

The  objective  of  this  study  is  to  investigate  numerically  a  compact  steam  methane  reforming  (SMR)  system  integrated  with  a  water-gas  shift 
(WGS)  reactor.  Separate  numerical  models  are  established  for  the  combustion  part,  SMR  and  WGS  reaction  bed.  The  concentration  of  species  at 
the  exits  of  the  SMR  and  WGS  bed,  and  the  temperatures  in  the  WGS  bed  are  in  good  agreement  with  the  measured  data.  Heat  transfer  to  the 
catalyst  beds  and  the  catalytic  reactions  in  the  SMR  and  WGS  catalyst  bed  are  investigated  as  a  function  of  the  operation  parameters.  The  conversion 
of  methane  at  the  exit  of  the  SMR  catalyst  bed  is  calculated  to  be  87%,  and  the  carbon  monoxide  concentration  at  the  outlet  of  the  WGS  bed  is 
estimated  to  be  0.45%.  The  effects  of  the  cooling  heat  flux  at  the  outside  wall  of  the  system  and  steam- to-carbon  (S/C)  ratio  are  also  examined.  As 
the  cooling  heat  flux  increases,  both  the  methane  conversion  and  carbon  monoxide  content  are  reduced  in  the  SMR  bed,  and  the  carbon  monoxide 
conversion  is  improved  in  the  WGS  bed.  Both  methane  conversion  and  carbon  dioxide  reduction  increase  with  increasing  steam-to-carbon  ratio. 
©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  fuel  cell  is  a  key  technology  for  realizing  a  hydrogen 
economy.  Various  types  of  fuel  cells  that  are  operated  with 
hydrogen  have  been  developed  to  date.  For  the  successful  com¬ 
mercialization  of  fuel  cell  a  stable  supply  of  low-cost  hydrogen 
is  required.  Although  many  technologies  have  been  developed, 
the  production  of  hydrogen  from  natural  gas  is  considered  to  be 
highly  promising. 

The  author’s  laboratories  have  conducted  studies  directed  at 
the  development  of  fuel  cells  and  hydrogen  production  with  an 
emphasis  on  polymer  electrolyte  membrane  fuel  cells  (PEM- 
FCs)  and  a  system  for  producing  hydrogen  from  natural  gas. 
The  PEMFC-based  cogeneration  systems  show  considerable 
promise  for  application  in  a  small-scale  distributed  power  gen¬ 
eration  in  households.  The  PEMFC  stack  is  integrated  with  a 
hydrogen-generation  system  to  provide  hot  water  and  electric¬ 
ity.  High  efficiency,  low  noise,  and  low  emission  of  air  pollutants 
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constitute  the  major  benefits  derived  from  the  use  of  a  PEMFC 
system  [1].  To  achieve  this  purpose,  a  small-scale  hydrogen  gen¬ 
eration  system  is  needed,  and  to  address  this  issue,  we  have 
developed  a  compact  steam  methane  reformer  using  natural  gas 
as  the  feedstock. 

Many  types  of  fuel-processing  system  for  producing  hydro¬ 
gen  that  is  sufficiently  pure  and  rapid  for  supplying  a  PEMFC  are 
currently  being  developed  [2-4] .  Although  the  steam-reforming 
process  requires  more  start-up  time  to  reach  a  steady-state  than 
an  autothermal  process,  it  can  yield  hydrogen  concentrations 
above  75%  in  the  dry  product  gas,  while  the  autothermal  pro¬ 
cess  produces  gas  containing  less  than  50%  hydrogen  because 
of  the  high  nitrogen  content  that  is  supplied  for  the  internal  com¬ 
bustion  of  natural  gas  and  air  [5].  Therefore,  a  natural  gas  fuel 
processor  in  which  a  steam-reforming  process  is  adopted  has 
been  widely  developed  as  one  of  the  key  components  of  the 
residential  PEMFC  system. 

The  present  study  focuses  on  developing  a  compact  fuel¬ 
processing  system  for  a  PEMFC  of  1.0  kW  capacity.  The  unit 
processes  of  a  gas  burner,  a  steam  methane  reforming  (SMR) 
reactor,  a  water-gas  shift  (WGS)  reactor  and  internal  heat 
exchangers  are  integrated  into  a  single  modular  unit.  Laboratory 
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Nomenclature 

Dp  Diameter  of  catalyst  particle  [m] 

E  Total  energy  [J] 

g  Gravitational  force  [ms-2] 

hi  Enthalpy  of  species  i  [Jkg-1] 

Ji  Diffusion  flux  of  species  i  [kgm-2  s-1] 
k  Thermal  conductivity  [wm-1  K-1], 

^eff  Effective  thermal  conductivity  [wm-1  K-1] 

kf  Thermal  conductivity  of  fluid  [w  m-1  K-1] 

ks  Thermal  conductivity  of  solid  [wm-1  K-1] 

k\,k2,  £3  Reaction  rate  constants  [see  Table  1] 

Ki  Constants  [see  Table  1] 

mi  Mass  fraction  of  species  i 

mp  Mass  fraction  of  product  species  p 

mR  Mass  fraction  of  reactant  species  R 

p  Pressure  [N  m-2] 

Pi  Partial  pressure  of  species  i  [bar] 

R  Universal  gas  constant  [8 . 3 1 4  J  mol- 1  K- 1  ] 

R,  Gas  reaction  rate  of  species  i  [mol  m-3  s-1] 
T  Temperature  [K] 

u  Velocity  [ms-1] 

Wi  Molecular  weight  of  species  i  [kg] 

Yi  Mole  fraction  of  species  i 

s  Dissipation  rate  of  turbulent  kinetic  energy 

0  Porosity 

rj i  Effectiveness  of  species  i 

k  Turbulent  kinetic  energy 

v'  Stoichiometric  coefficient  for  reactant 

v"  Stoichiometric  coefficient  for  product 

p  Density  [kgm-3] 

r  Stress  tensor 


Table  1 

Parameters  for  kinetic  model  for  methane  steam  reforming 

Pre-exponential  Activation 

factor  (Aj)  energy  (Ep 

[dimension  of  k]  [kJ  mol  - 1  ] 


Reaction  rate  constants  for  Arrhenius 
expression:  ki  =A;exp(— Zq/RT ) 

h  [kmol  bar0  5  kg  cat-1  h-1]  4.225  x  1015  240.1 

k2  [kmol  bar-1  kg  cat-1  h-1]  1.955  x  106  67.13 

k3  [(kml  bar0  5  kg  cat-1  h-1)]  1.020  x  1015  243.9 


Pre-exponential 
factor  (Bp 
[dimension  of  K] 

Enthalpy  change  of 
adsorption  (AHp 
[kJ  mol-1] 

Constants  for  Van’t  Hoff  equation: 

Ki  =  BiQx  p(- 

-AHJRT) 

^co  [bar-1] 

8.23  x  10“5 

-70.61 

Ku2  [bar-1] 

6.12  x  10"9 

-82.90 

^CH4[bar-1; 

|  6.65  x  10“4 

-38.28 

Kn2o 

1.77  x  105 

88.68 

testing  is  performed  to  investigate  the  characteristics  of  the  unit 
and  to  improve  its  performance  further. 

A  numerical  simulation  is  a  very  useful  tool  for  developing  an 
efficient  SMR  system.  Using  numerical  simulation,  the  effects  of 
design  and  operating  parameters  can  be  investigated  at  a  low  cost 
and  the  results  can  provide  basic  data  for  designing  and  experi¬ 
mentally  testing  such  systems.  In  this  study,  the  SMR  system  is 
numerically  investigated  with  respect  to  heat  transfer  between 
the  heat  source  and  the  reacting  catalyst  bed,  the  reactions  in 
the  SMR  and  WGS  catalyst  beds,  the  effects  of  operating  con¬ 
ditions,  etc.  To  accomplish  this,  a  numerical  model  applicable 
to  the  reforming  system  is  established  and  the  models  are  val¬ 
idated  by  comparing  numerical  results  with  experimental  data. 
After  validation,  the  reforming  system  is  analyzed  with  respect 
to  the  various  effects  of  operation  and  design  parameters  on  heat 
transfer  performance  and  catalytic  reaction. 

2.  Numerical  models 

The  numerical  simulation  for  the  SMR  system  is  carried 
out  by  separating  it  into  three  different  parts,  based  on  its 
structure  (see  Fig.  1).  One  part  corresponds  to  the  heat  sup¬ 
ply  which  is  comprised  of  a  gas  burner,  a  combustion  chamber 
and  a  channel  for  the  combustion  gas.  The  second  part  is  the 
SMR  reaction  region  of  the  feedstock,  which  contains  reform¬ 
ing  catalysts.  The  third  represents  a  WGS  reactor  that  comprises 
a  high-temperature  shift  (HTS)  reactor  (300-400  °C)  and  a 
low-temperature  shift  (LTS)  reactor  (200-300  °C).  A  separate 
numerical  model  was  established  for  each  part.  The  following 
describes  each  of  the  numerical  models  in  detail. 

2.7.  Model  of  combustion  and  heat  transfer 

The  numerical  models  described  here  include  the  modelling 
of  combustion  in  the  combustion  chamber  and  heat  transfer  asso¬ 
ciated  with  the  combustion  gas.  The  computational  domain  was 
set  up  as  an  axisymmetric  system  that  reflects  the  fact  that  the 
reforming  system  has  a  cylindrical  shape.  The  governing  equa¬ 
tions  for  heat  generation  and  heat  transfer  can  be  written  as 
follows. 

Continuity  equation: 

V  •  (pu)  =  0  (1) 

LTS  Water  Tube 

Porous  Media 


NG+H20 


NG+Air 


Fig.  1.  Schematic  of  compact  SMR  system  incorporated  with  WGS  reactor. 
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Momentum  equation: 

V  •  ( puu )  =  -Vp  +  Vr  +  p  g 


combustion  gas.  The  WSGGM  is  regarded  to  be  a  quite  reason¬ 
able  model  in  combustion  applications  and  detailed  explanations 
(2)  can  be  found  in  the  literature  [10,11]. 


Energy  equation: 


V  •  (u(pE  +  P))  =  V  •  (&effVr  +  ^2  h0i  +  (“  •  reff)) 


+  hiRi  +  V  •  5rad 


(3) 


Species  equation: 

V  •  (p  u  mi)  =  -  V  •  },  +  (4) 

State  equation: 

e  =  “RT'Ew,  <5) 

In  the  above  governing  equations,  turbulent  flow  is  modelled 
using  the  RNG-based  k-s  two-equation  model.  Details  of  the 
RNG  model  can  be  found  in  the  literature  [6].  The  RNG  k-s 
model  is  known  to  be  more  accurate  and  reliable  than  the  stan¬ 
dard  k-s  turbulent  models  for  a  wide  class  of  flows.  The  com¬ 
bustion  in  the  gas  burner  is  modelled  using  the  finite-chemistry 
reaction  model  and  the  eddy-dissipation  model  (EDM)  [7] .  Most 
fuels  that  are  used  in  gas  combustion  are  fast  burning,  and  the 
overall  rate  of  such  reactions  can  be  assumed  to  be  controlled 
by  turbulent  mixing.  In  such  a  pre-mixed  flame,  as  in  the  present 
study,  the  turbulence  slowly  mixes  cold  reactants  and  hot  prod¬ 
ucts  into  the  reaction  zone,  in  which  the  reaction  starts  to  occur 
rapidly.  In  such  a  situation,  the  combustion  is  said  to  be  mixing- 
limited,  and  the  complex  chemical  kinetic  rate  can  be  safely 
neglected.  Thus,  the  use  of  the  eddy-dissipation  model  pro¬ 
vides  the  advantage  of  a  computationally  low  cost  approach. 
This  model  calculates  the  net  production  rate  with  the  smaller 
of  the  two  following  equations: 

.  s  (  wr  \ 

Ri  =  VfWiAp- mm  (  — —  )  (6) 

k  \vRWR ) 


2.2.  Model  of  SMR  reaction 


The  SMR  catalyst  bed  is  modelled  as  a  porous  media,  which 
reflects  the  fact  that  a  reactor  is  packed  with  spherical  cata¬ 
lysts  of  diameters  less  than  5  mm.  Homogeneity  is  assumed  in 
calculating  the  energy  equation  for  the  SMR  catalyst  bed.  The 
homogeneous  model  assumes  that  the  temperatures  of  the  cata¬ 
lysts  are  the  same  as  those  of  the  fluid  around  the  catalysts.  The 
governing  equations  for  the  SMR  catalyst  bed  can  be  written  as 
follows: 

Continuity  equation: 

V  •  (pu)  =  0  (8) 


Momentum  equation: 


V  •  (puu)  =  — V p  +  Vr  ■ 


150(1  -0)2 


•  / iu 


Dj  03 

Energy  equation: 

V  •  (u(pE  +  p))  =  V  •  (keffVT  +  ^2  hiJi  +  (u  ■  reff)) 
+  ^2  hiViK 


(9) 


(10) 


where  keR  =  cp  K(  +  ([+  <P)ks 
Species  equation: 

V  •  (p  u  mi)  =  -V  •  Jt  +  rji R,  (11) 


State  equation: 

p  =  i>rtY.w,  <12) 

The  above  energy  and  species  equations  require  the  calculation 
of  the  SMR  reaction  in  the  catalyst  bed.  The  SMR  reaction  can 
be  expressed  by  the  following  equations: 


Ri  =  v-WiABp 


K^v'/Wi 


(7)  CH4  +  H20  CO  +  3H2  Zlffi98  =  206  kJ  moF 1 

(13) 


where  A  and  B  are  empirical  constants  with  values  of  4.0  and 
0.5,  respectively. 

The  flame  produced  by  gas  combustion  maintains  a  very  high 
temperature,  over  1000  °C,  where  radiation  becomes  a  signifi¬ 
cant  factor  in  heat  transfer.  The  combustion  gas  participates  in 
the  radiation  process,  that  is,  it  absorbs,  emits  and  scatters  radia¬ 
tion  energy.  Such  a  situation  requires  modelling  of  the  radiative 
heat  transfer  for  the  combustion  gas.  The  present  study  employs 
the  discrete  transfer  radiation  model  (DTRM)  [8,9]  for  the  radia¬ 
tive  heat  transfer  of  the  combustion  gas. 

In  calculating  the  radiation  of  combustion  gas,  the  selection 
of  the  appropriate  radiative  thermal  properties,  one  of  which  is 
the  absorption  coefficient  of  the  combustion  gas,  is  an  important 
issue.  The  present  modelling  uses  the  weighted-sum-of-grey- 
gases-model  (WSGGM)  for  the  absorption  coefficient  of  the 


CO  +  H20  -o-  C02  +  H2  Z\H298  =  -4 1 . 1  kJ  moP 1 

(14) 


CH4  +  2H20  C02  +  4H2  Z\H298  =  164.9  kJ  moP 1 

(15) 


The  kinetics  of  the  SMR  reaction  have  been  studied  by  several 
researchers  [12-14].  Xu  and Froment’s  results  [12]  are  employed 
in  the  present  study.  The  intrinsic  rate  equations  proposed  by 
these  authors  are  as  follows: 


n 


k\  (pcnAPn2o  ~  Ph2Pco/K\) 

pU  (DEN? 


(16) 
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h.  (pcoPHiO  -  Pn2 Pah / Ko ) 
PH2  (DEN)2 

h  (Pch4Ph2o  -  pfpcOi/K-1) 

pho  (DEN)2 


tested  individually  before  coupling.  The  results  of  the  numeri¬ 
cal  calculation  are  compared  with  experimental  data  in  order  to 
validate  the  models. 


(18) 


3.  Results  and  discussion 


with  the  denominator: 


DEN  =  1  +  K  coPco  +  Kh2Ph2  +  ^ch4/?ch4  + 


Kh2oPu2o 

Pn2 

(19) 


The  activation  energies  and  pre-exponential  factors  for  the 
reaction  rate  constants  k\,  and  £3,  the  adsorption  enthalpy 
change  and  pre-exponential  factors  for  the  adsorption  constants 
^co ,  KU2 ,  ^ch4  and  ^h2o  are  listed  in  Table  1 . 

For  calculating  the  reforming  reaction  in  the  SMR  catalyst 
bed,  the  effectiveness  factor,  77,  in  the  above  Eqs.  (10)  and  ( 1 1 )  is  a 
significant  parameter.  The  present  model  uses  the  assumption  of 
homogeneity  between  fluid  and  catalysts,  and  the  effectiveness 
factor  is  not  calculated  separately,  but  is  selected  as  0.05  on  the 
basis  of  effectiveness  analyses  carried  out  by  others  [15,16]. 


2.3.  Model  ofWGS  reaction 


The  WGS  catalyst  bed  is  a  porous  media,  the  same  as  the 
previous  SMR  reforming  catalyst  bed.  Therefore,  the  govern¬ 
ing  equations  for  the  WGS  catalyst  bed  are  the  same  as  those 
for  the  SMR  catalyst  bed  except  for  the  kinetics  of  the  WGS 
reaction.  In  general,  the  WGS  catalyst  bed  consists  of  two  oper¬ 
ation  zones,  an  HTS  bed  operated  at  300-400  °C  and  an  LTS 
bed  operated  at  200-300  °C.  Fe304/Cr203  and  Cu/Zn0/Al203 
are  typically  used  as  catalysts  for  the  high-and  low-temperature 
bed,  respectively.  The  present  calculation  uses  the  same  kinetics 
for  the  HTS  and  LTS  beds  for  ease  of  calculation. 

The  WGS  reaction  of  the  CO  component  can  be  expressed 
by  the  following  equation. 

CO  +  H20  CO2  +  H2  ziH°98  =  —41.1  kJ  mol-1 

(20) 


The  kinetics  have  been  studied  by  several  researchers  [17-21]. 
Moe’s  results  [20]  are  employed  in  the  present  study,  the  rate 
equation  is  as  follows: 


y 4  =  k4  pco  Pn2o 


Pu2Pco2  \ 
K2PC0PU20  ) 


(21) 


<  (  1 855.5  \  1 

k\  =  1.85  x  10  5exp  (  12.88 - — —  j  [molg  1  min  l] 

(22) 

To  solve  the  complete  set  of  governing  equations  established 
for  gas  combustion,  the  SMR  reaction  and  the  WGS  reaction, 
the  FLUENT  software  program  [22]  is  used.  Since  three  parts 
of  gas  combustion,  SMR  and  WGS  reactions,  are  coupled  with 
each  other,  their  governing  equations  are  solved  simultaneously. 
To  ensure  the  accuracy  of  the  calculation,  however,  each  part  is 


3.1.  Characteristics  of  combustor,  SMR  and  WGS  reactor 

The  SMR  reaction  is  a  highly  intensive  endothermic  process. 
The  heat  required  for  the  reforming  reaction  is  supplied  by  the 
combustion  gas  via  a  gas  burner.  A  schematic  of  the  SMR  system 
integrated  with  a  WGS  reactor  is  shown  in  Fig.  1.  For  effective 
heat  transfer  from  the  combustion  gas  to  the  SMR  catalyst  bed, 
a  narrow  flow  channel  is  placed  next  to  the  SMR  catalyst  bed. 
The  combustion  gas  supplies  the  main  heat  energy  to  the  SMR 
catalyst  bed  during  its  passage  through  the  channel.  A  mixture 
of  natural  gas  (CH4)  and  steam,  as  a  feedstock,  is  fed  to  the 
SMR  catalyst  bed,  which  absorbs  the  heat  energy  from  the  com¬ 
bustion  gas.  Before  the  feedstock  reaches  the  SMR  catalyst,  it 
passes  through  porous  media,  which  fills  the  channel  to  produce 
a  uniform  flow  of  feedstock  and  also  to  enhance  the  heat  trans¬ 
fer  from  the  combustion  gas  to  the  feedstock.  The  feedstock  in 
the  SMR  catalyst  is  converted  into  reformed  gas,  which  then 
flows  into  the  WGS  catalyst  bed  (HTS  +  LTS),  located  in  the 
outermost  channel  of  the  system.  The  WGS  bed  reduces  the  car¬ 
bon  monoxide  content  of  the  reformed  gas  to  less  than  1.0%. 
The  WGS  reaction  generates  a  slight  amount  of  heat  and  must 
be  cooled.  To  accomplish  this,  cooling  water  tubes  are  installed 
around  the  WGS  bed.  The  cooling  water  in  the  tubes  is  converted 
into  steam  by  the  heat  transfer  and  is  used  as  the  feedstock  of 
the  SMR. 

To  check  the  validity  of  the  numerical  simulation,  the  cal¬ 
culated  results  were  compared  with  experimental  data.  Fig.  2 
illustrates  a  comparison  between  the  experimental  and  calcu¬ 
lated  on  temperatures  of  the  WGS  catalyst  bed.  The  calculated 
temperatures  show  slightly  lower  values  in  the  inlet  of  the  WGS 


Fig.  2.  Comparison  of  calculated  and  measured  temperatures  along  centerline 
of  WGS  bed. 
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Fig.  3.  Comparison  of  calculation  and  measured  data  on  species  at  (a)  exit  of 
SMR  catalyst  bed  and  (b)  WGS  catalyst  bed. 


bed  and  are  slightly  higher  in  the  exit  than  the  measured  data.  A 
comparison  of  the  products  in  the  SMR  and  WGS  catalyst  bed 
is  shown  in  Fig.  3.  The  calculated  values  for  CH4  and  H2  are 
very  close  to  the  experimental  data,  while  the  carbon  monox¬ 
ide  and  CO2  contents  deviate  slightly  from  the  measured  data. 
Such  a  difference  may  be  caused  by  an  inaccurate  calculation  of 
reaction  rate  or  any  other  reason.  Even  though  some  discrepancy 
exists  in  the  calculation  results,  the  numerical  simulation  would 
be  expected  to  provide  a  useful  tool  for  analyzing  the  effects  of 
various  design  and  operation  parameters. 

Heat  transfer  from  the  combustion  gas  to  the  SMR  and  WGS 
catalyst  bed  was  investigated.  Fig.  4  illustrates  the  temperature 
distribution  in  the  combustor  and  the  catalyst  beds  of  SMR  and 
WGS.  Within  the  combustor,  the  highest  temperature  is  along  the 
central  line  of  the  combustor,  which  could  alleviate  the  risk  of  hot 
spots  on  the  wall  of  the  SMR  catalyst  bed.  Nevertheless,  the  top 
end  of  the  combustor  reaches  a  very  high  temperature  and  should 
be  made  of  sufficiently  strong  steel  plate.  The  combustion  gas 
passes  through  the  channel  between  the  combustor  and  the  SMR 


1 


1  939*03 
1  860*03 
1  8le*03 
1  76e*03 
'  720*02 
1  67e*03 
1  629*03 
1  57e*03 
1.520*03 
1  47e*03 
1  42o*03 
1  370*03 
1  33e*03 
1  28e*03 
1  230*03 


1  130*03 
1  080*03 


I 


9  830*02 
9  35e*02 
8.86e*02 
8  370*02 


7880*02 
7  390*02 
6  9-0*02 


642e*02 
5 93©* 02 
5  44e*02 
4.95e*02 
4  460*02 
3.980*02 
3  490*02 
3.000*02 


(a) 


I 

I 


1  93«*03 
1  860*03 
1  81e*03 
1  750*03 
1  72e*03 
1  67e*03 
1  620*03 
1  570*03 
1  520*03 
1  470*03 
1  420*03 
1  370*03 
1  330*03 
1  280*03 
1  23«*03 
1  180*03 
1  130*03 
1  080*03 
1  030*03 
9  8  30  *  02 
9  350*02 
8  860*02 
8  370*02 


78 


>*02 


7  39e*02 
6  910*02 
6  4  20  *  02 
5  930*02 
5  440*02 
4  950*02 
4  4  60*02 
3  980*02 
3  490*02 
3  00e*02 


(b) 


Fig.  4.  (a)  Temperature  and  (b)  path  line  distribution  in  SMR  system. 


catalyst,  and  the  heat  is  transferred  from  the  combustion  gas 
into  the  SMR  bed  for  methane  steam  reforming,  an  endothermic 
reaction. 

Fig.  5  describes  the  temperature  profile  of  the  combustion 
gas,  in  the  channels  of  the  SMR  and  WGS  bed.  In  the  channel 
for  the  combustion  gas,  the  temperature  drops  rapidly  from  about 
1 100  °C  to  below  320  °C  along  the  flow  of  combustion  gas.  This 
is  due  to  heat  transfer  from  the  combustion  gas  to  the  SMR 
catalyst  bed.  The  temperature  in  the  SMR  catalyst  bed  next  to 
the  channel  for  the  combustion  gas  increases  linearly  from  about 


Fig.  5.  Temperature  profiles  along  combustion  gas  channel,  SMR  and  WGS 
catalyst  bed. 
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100  °C  to  around  640  °C.  The  temperature  in  the  WGS  channel 
decreases  steadily  from  616  °C  at  the  inlet  to  190  °C  at  the  outlet. 
Such  a  decrease  in  the  temperature  of  the  WGS  bed  is  caused 
by  the  cooling  water  tubes  installed  around  the  WGS  channel. 
The  computed  data  show  that  temperatures  in  the  channel  for 
the  combustion  gas  are  by  100-250  °C  higher  than  those  in  the 
SMR  bed,  but  the  temperatures  in  the  SMR  bed  are  nearly  the 
same  as  those  in  the  WGS  bed.  The  reason  for  this  is  that  the 
SMR  bed  and  the  WGS  bed  are  filled  with  catalysts,  which  act  as 
media  to  augment  heat  transfer.  Similarly  if  any  porous  medium 
is  installed  in  the  channel  for  the  combustion  gas,  it  will  improve 
the  heat  transfer  to  the  SMR  bed  from  the  combustion  gas. 

Fig.  6  shows  the  heat  flux  on  to  the  wall  of  the  SMR  catalyst 
bed.  The  heat  flux  is  in  the  range  5.7-16.0  kW  m-2  in  the  region 
of  the  porous  media  and  16.0-55.2  kW  m-2  in  the  region  of  the 
SMR  catalyst.  The  amount  of  heat  transferred  from  the  com¬ 
bustion  gas  is  calculated  to  be  3.08  kW  in  the  porous  media  and 
10.46  kW  in  the  SMR  catalyst.  The  total  heat  transferred  from 
the  combustion  gas  to  the  SMR  channel  is  13.54 kW,  which  is 
estimated  to  be  76.1%  of  the  total  heat  energy  supplied  from 
the  gas  burner.  The  residual  heat  contained  in  the  flue  gas  can 
be  recovered  through  additional  heat  exchangers  (not  shown  in 
Fig.  1.)  Achieving  a  high  thermal  efficiency  is  an  important  tar¬ 
get  in  developing  a  new  compact  reformer  system.  The  higher 
heat  recovery  from  the  combustion  gas  leads  to  a  more  efficient 
reformer  system  in  terms  of  thermal  energy.  The  analysis  sug¬ 
gests  that  it  would  be  desirable  to  enhance  further  heat  transfer 
from  the  combustion  gas  to  the  SMR  channel. 

The  concentration  profiles  of  species  within  the  SMR  catalyst 
bed  are  depicted  in  Fig.  7(a).  The  CH4  concentration  decreases 
rapidly  from  the  inlet  of  the  SMR  catalyst  and  reduces  to  its 
lowest  level  (2.8%)  at  the  outlet  of  the  catalyst  bed.  The  H2  con¬ 
centration  rapidly  increases  from  the  inlet  and  reaches  a  peak 
near  the  outlet.  The  profile  for  H2  concentration  is  inversely 
analogous  to  the  CH4  concentration  profile.  By  contrast,  pattern 
for  the  CO  concentration  profile  is  slightly  different  from  that 
for  H2.  The  CO  concentration  in  the  front  part  of  the  SMR  bed 
increases  more  slowly  than  in  the  rear  part  of  the  SMR  bed.  On 
the  other  hand,  Fig.  7(b)  represents  the  reaction  rate  for  each  of 


the  reaction  process,  rl,  r2  and  r3  in  Eqs.  (13)— (15).  The  data 
reveal  that  the  rl  reaction  is  the  dominant  process  among  the 
three  reaction  processes.  Its  rate  increases  from  the  inlet  of  the 
SMR  bed  and  reaches  the  maximum  after  the  middle  of  the  bed, 
then  decreases  rapidly.  The  r2  reaction  shows  a  backward  reac¬ 
tion  near  the  exit  of  the  bed,  depicted  as  a  negative  reaction  rate. 
The  r2  reaction  is  a  mild  exothermic  reaction  and  the  equilibrium 
constant  decreases  with  increasing  temperature.  As  described 
above,  the  temperature  in  the  SMR  bed  increases  along  the  flow 
of  reactants.  The  r2  reaction  reaches  a  peak  around  the  middle  of 
the  catalyst  bed,  then  decreases  steadily,  and  eventually  reaches 
a  negative  value. 

The  reformed  gas  exiting  from  the  SMR  catalyst  bed  is  then 
fed  to  the  WGS  catalyst  bed,  where  the  CO  is  reduced  via  the 
WGS  reaction  (Eq.  (20)).  Fig.  8  represents  the  species  profiles 
and  the  WGS  reaction  rate.  The  CO  decreases  rapidly  from  8.7% 
at  the  inlet  to  0.5%  at  the  outlet,  while  the  CO2  increases  from 
12.5  to  19.1%  through  the  WGS  bed.  The  CH4  concentration 
shows  little  change  through  the  catalyst  bed.  Fig.  8(b)  shows 
a  plot  of  the  reaction  rate  of  Eq.  (20)  in  the  WGS  bed.  The 
reaction  rate  decreases  drastically  near  the  inlet  and  then  slowly 
decreases  toward  the  outlet.  The  temperature  in  the  WGS  bed 
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Fig.  8.  (a)  Profiles  of  species  and  (b)  reaction  rate  in  WGS  catalyst  bed. 


changes  from  489  °C  at  the  inlet  to  190  °C  at  the  outlet,  as  the 
result  of  cooling  at  the  outside  of  the  WGS  bed.  The  heat  flux  is 
set  to  7000  W  vcT1  at  the  outside  wall  of  the  WGS  bed. 

3.2.  Effects  of  operation  parameters 

For  a  specified  SMR  system,  several  operation  parameters 
can  significantly  influence  the  performance  of  the  system.  One 
of  them  is  thought  to  be  a  cooling  heat  flux  at  the  outermost 
wall  of  the  system.  The  cooling  lowers  the  temperature  of  the 
WGS  catalyst  bed  so  that  it  operates  at  a  low  temperature  for  the 
efficient  conversion  of  CO.  As  described  above,  the  HTS  is  oper¬ 
ated  at  temperatures  of  300-400  °C  and  the  LTS  at  temperatures 
of  200-300  °C.  The  amount  of  the  cooling  heat  flux  would  be 
expected  to  affect  directly  not  only  the  performance  of  HTS  and 
LTS  but  also  the  performance  of  the  SMR  bed.  Fig.  9  depicts  the 
effects  of  a  cooling  heat  flux  on  temperatures  in  the  WGS  bed. 
These  decrease  with  increasing  cooling  heat  flux.  The  temper¬ 
atures  near  the  exit  of  the  WGS  bed  are  more  heavily  affected 
by  the  cooling  heat  flux  than  those  near  its  inlet.  Fig.  10  shows 
the  effects  of  the  cooling  heat  flux  on  the  reforming  reaction  in 
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Fig.  9.  Effects  of  cooling  heat  flux  at  outermost  wall  of  system  on  temperatures 
in  WGS  bed. 
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Fig.  10.  (a)  Species  at  the  SMR  outlet  and  (b)  WGS  outlet  as  function  of  cooling 
heat  flux  at  outermost  wall  of  system. 
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the  SMR  and  on  CO  conversion  in  the  WGS  bed.  As  the  cool¬ 
ing  heat  flux  increases,  methane  conversion  and  the  CO  content 
in  the  SMR  bed  both  decrease.  The  reason  for  the  decrease  in 
methane  conversion  is  that  the  higher  cooling  heat  flux  lowers 
the  temperature  of  the  SMR  bed,  so  that  the  reactions  of  Eqs. 
(13)— (15)  move  to  an  equilibrium  corresponding  to  the  lower 
temperature  and  thus  result  in  a  lower  conversion  of  methane. 
This  suggests  that  the  cooling  heat  flux  should  be  controlled 
appropriately  so  as  not  to  drop  the  conversion  of  methane  to  less 
than  the  target  level.  As  an  example,  the  target  for  CH4  con¬ 
version  in  the  present  development  is  over  85%,  which  requires 
the  cooling  heat  flux  to  be  maintained  at  less  than  7000  W  m-2. 
An  investigation  concerning  the  CO  conversion  in  the  WGS  bed 
reveals  that  it  increases  as  the  cooling  heat  flux  is  increased  (see 
Fig.  10(b)).  The  CO  at  the  outlet  of  the  WGS  bed  is  also  reduced 
with  increase  in  the  cooling  heat  flux.  This  analysis  suggests  that 
the  cooling  heat  flux  should  be  as  high  as  possible  to  reduce  the 
CO  content  at  the  exit  of  the  WGS  bed.  The  maximum  cooling 
heat  flux  is  limited  due  to  the  target  for  methane  conversion.  In 


Fig.  11.  Effects  of  S/C  ratio  (a)  reforming  reaction  in  SMR  and  (b)  CO  conver¬ 
sion  in  WGS  bed. 


the  case  of  the  present  design,  it  is  recommended  that  the  cooling 
heat  flux  be  increased  up  to  7000  W  m-2  to  obtain  the  maximum 
conversion  of  CH4  in  the  SMR  bed  and  the  maximum  reduction 
of  CO  in  the  WGS  bed. 

Another  important  parameter  of  the  reforming  system  is  the 
steam-to-carbon  (S/C)  ratio  of  the  feedstock.  The  steam  con¬ 
tent  in  the  feedstock  is  very  important.  The  theoretical  value  for 
the  S/C  ratio  is  1.0,  as  shown  in  Eq.  (13),  however,  in  practice 
coking  of  carbon  occurs  in  the  reforming  catalyst  of  the  SMR 
bed.  Accordingly  the  S/C  ratio  is  maintained  over  1 .0  to  avoid 
carbon  coking  and  also  to  increase  the  conversion  of  methane. 
Thus,  the  effects  of  the  S/C  ratio  on  the  reforming  reaction  in  the 
SMR  bed  and  on  CO  conversion  in  the  WGS  bed  were  inves¬ 
tigated;  the  results  are  plotted  in  Fig.  11.  The  results  for  the 
SMR  bed  show  that  methane  conversion  increases  as  the  S/C  is 
increased,  leading  to  an  increase  in  the  content  of  H2  and  CO2 
and  a  decrease  in  CO  content.  The  S/C  ratio  also  significantly 
affects  CO  conversion  in  the  WGS  bed,  as  shown  in  Fig.  11(b). 
As  the  S/C  ratio  increases,  the  CO  conversion  is  enhanced.  The 
calculation  also  shows  that  the  degree  of  enhancement  in  CO 
conversion  is  stronger  than  that  for  methane  conversion.  The 
analysis  implies  that  the  S/C  ratio  should  be  maintained  as  high 
as  possible.  On  the  other  hand,  maximum  supply  of  steam  is 
limited  within  the  allowable  cost  of  steam  generation.  A  higher 
S/C  ratio  leads  to  a  lower  thermal  efficiency  of  the  SMR  system. 
Therefore,  the  S/C  should  be  properly  determined  considering 
the  relation  between  the  increase  in  the  conversion  of  methane 
and  CO  and  the  disadvantage  in  thermal  efficiency. 

4.  Conclusions 

To  investigate  a  compact  SMR  system  integrated  with  a  WGS 
reactor,  a  numerical  simulation  has  been  conducted.  For  the 
numerical  calculation,  numerical  models  were  first  established. 
The  present  SMR  system  consists  of  three  parts:  a  combustor, 
an  SMR  catalyst  bed,  and  a  WGS  catalyst  bed.  Thus,  the  numer¬ 
ical  model  is  separately  set  up  for  each  part.  One  is  a  model  for 
combustion  and  heat  transfer,  and  the  other  two  are  models  for 
the  SMR  and  WGS  reactions.  The  SMR  and  WGS  catalyst  beds 
are  modelled  as  a  porous  media  and  are  assumed  to  be  homoge¬ 
neous  in  calculating  the  energy  equation.  After  establishing  the 
models,  the  calculations  are  validated  by  comparing  the  calcu¬ 
lation  results  with  experimental  data.  The  species  concentration 
at  the  exit  of  the  SMR  and  WGS  bed  and  the  temperatures 
in  the  WGS  bed  are  in  good  agreement  with  the  measured 
data. 

After  establishing  the  numerical  model,  the  SMR  system  is 
investigated  with  respect  to  heat  transfer  to  the  catalyst  beds, 
and  catalytic  reactions  in  the  SMR  and  WGS  catalyst  bed.  The 
amount  of  heat  transferred  to  the  SMR  bed  from  the  gas  burner 
is  estimated  to  be  76.1%  of  the  total  heat  energy  supplied  from 
the  combustion  gas.  The  calculation  suggests  that  the  heat  trans¬ 
ferred  to  the  catalyst  bed  should  be  further  increased  in  order  to 
enhance  the  thermal  efficiency  of  the  SMR  system.  The  conver¬ 
sion  of  CH4  at  the  exit  of  the  SMR  catalyst  bed  is  calculated  to 
be  87%,  and  CO  at  the  outlet  of  the  WGS  bed  is  estimated  to  be 
0.45%.  It  is  also  shown  that  for  the  r2  reaction  (Eq.  (14)),  the 
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backward  reaction  takes  place  in  the  high-temperature  region  of 
the  SMR  bed. 

The  effects  of  operation  parameters,  such  as  cooling  heat  flux 
at  the  outside  wall  of  the  system  and  steam-to-carbon  (S/C)  ratio, 
have  also  been  investigated.  As  the  cooling  heat  flux  increases, 
both  CH4  conversion  and  CO  content  are  reduced  in  the  SMR 
bed,  and  CO  conversion  is  improved  in  the  WGS  bed.  The 
analysis  implies  that  the  cooling  heat  flux  can  be  increased  to 
7000  Wm-2  while  CH4  conversion  is  maintained  at  over  85%. 
On  the  other  hand,  the  results  show  that  the  S/C  ratio  has  a 
significant  influence  on  CH4  conversion  in  the  SMR  and  CO 
reduction  in  the  WGS.  As  the  S/C  ratio  is  increased,  CH4  con¬ 
version  and  CO  reduction  both  increase.  The  S/C  ratio  should 
be  maintained  at  as  high  a  value  as  possible,  but  the  maximum 
supply  of  steam  is  limited  within  the  allowable  cost  of  steam 
generation. 
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